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ABSTRACT 


The cross-spectra of temperature and ozone mass mixing ratio at 42 km 
{2 mb) and 28 km (20 mb) have been determined for austral spring (1971) 
and summer (1971-2) over Christchurch, New Zealand (44 S, 172 E). Tlie 
sources of data are the SCR and BUV experiments on Nimbus 4. The 
observed covariances are compared with a model in which the temperature 
and ozone perturbations are forced by an upward-propagating planetary 
wave. The agreement bebveen the observations and the model is reason- 
able. It is suggested that this cross- spectral method permits an estimate 
of the meridional gradient of ozone mass mixing ratio from measurements 
of the vertical profile of ozone mass mixing ratio at one location, supported 
by temperature profiles from at least two locations (to determine the 
meridional temperature gradient). 
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1. Introduction 


The availability of global measurements of stratospheric temperature and 
ozone from satellites is increasing our knowledge of the horizontal, vertical and 
temporal distribution of ozone. However, balloon-borne instrumentation is still 
a major source of information, particularly for the period before extensive 
satellite observations became available. It is, therefore, interesting to consider 
how much information on atmospheric motion can be deduced from the various 
time series of balloon observations. 

The association between ozone and atmospheric temperature is commonly 
described in terms of eddy transport coefficients, with the implication that 
dynamic processes with time scales less than several months, are not measured 
directly and are only represented by their integrated effects. The balloon 
observations are usually made daily so that motion with periods of two or more 
days can be satisfactorily measured. It is, therefore, possible to study the ozone 
transport on quite short time scales and, hence, obtain an insight into the dy- 
namical processes contributing to the eddy transport coefficients. 

The temperature and the ozone concentration at a fixed point are influenced 
considerably by advection of horizontal gradients so that information on these 
gradients is essential. Altliough the relevant temperature gradients can usually 
be obtained from a group of radiosonde stations, the ozone mixing ratio gradient 
is seldom available. 
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This investigation simulates a set of balloon observations by using the 
satellite measurements of temperature and ozone mixing ratio at a fixed geo- 
graphic location (Christchurch, New Zealand, 44 S, 172 E). The cross spectruni 
of the daily temperature and ozone mixing ratio measurements are calculated 
for austral spring (1971) and summer (1971-2). Suitable ozone data for other 
seasons were not available at the time of this study. The cross spectrum is 
then compared with that calculated from a simple model based on specified 
planetary wave modes. The observed seasonal mean values were used to define 
the atmospheric properties for the model. An attempt was then made to decide 
whether the meridional mixing ratio gradient could be deduced from the phase 
spectrum of the temperature-ozone cross -spectra. 

2. The Observations 

Temperature. Temperature information was talcen from the Nimbus 4 SCR 
(Barnett et al. , 1972) radiance tapes available from the Department of Atmo- 
spheric Physics, University of Oxford. Tlie black-body equivalent of the radiance 
is used, i.e. the temperature of a black body which would emit the same radiance 
at the same wavelength. The Nimbus 4 SCR weighting functions for the two top 
channels have maxima near 42 km (2 mb) and 28 km (20 mb). Tlie equivalent 
temperatures are weighted means over altitude regions about 20 km thick. 

The Nimbus 4 satellite is sun-synchronous so that the observations are 
made near local noon and local midnight. The data used in the analysis was 
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obtained from the "gridded” tapes where all the radiances for one day have been 
reduced to a latitude-longitude grid, the grid points being 4" apart in latitude 
and 10“ apart in longitude. The equivalent temperatures were obtained from 
daily radiances at the grid point (44 S, 172 E) closest to Christchurch, New 
Zealand. The mean meridional temperatui’e gradient was calculated from the 
difference between the grid points (40 S, 172 E) and (48 S, 172 E). 

Ozone. The ozone data was obtained from the Nimbus 4 BUV (Heath et al., 
1973) results available at NASA/Goddard Space Flight Center. Vertical profiles 
of ozone partial pressure were given at successive locations along orbits pass- 
ing near Christchurch. This BUV data is to be considered as a preliminary 
estimate until some inconsistencies in the data have been resolved (Krueger 
et al., 1973). 

As the SCH weighting functions peak at Ian (2 mb) and 28 Ian (20 mb), the 
ozone partial pressures for these heights were selected. The .mean value and 
the meridional gradient were obtained by fitting a least squares plane to the data 
fx'om a minimum of two orbits, one on either side of Christchurch. On the few 
occasions when only one orbit’s data was available for a given day, the mean 
value and the gradient were derived from a least squares line. In this latter . 
case there is the possibility of contamination due to zonal temperature gradients. 
Only data points taken within 3300 km of Christchurch were used, Tlie ozone 
mass mixing ratio was derived as the ratio of ozone partial pressure to the total 
air pressure, multiplied by the ratio of the respective molecular weights (1.66). 
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Vertical Resolution. The vertical resolution of the satellite experiments 


is limited by the width of the weighting functions. The half-weight width of the 
weighting functions for the upper two Nimbus 4 SCR channels have been estimated 
from the published (Barnett et al., 1972) weighting function graphs to be 20 km. 
Similarly the BUV half-weight widths are 17-18 km for the upper five channels 
(Heath et al., 1973), These weighting functions will only slightly attenuate the 
response to a wave with a vertical wavelength whose half-wavelength is close to 
the half-weight width, i.e. a wave whose wavelength is 35-40 km. Waves of 
shorter vertical wavelengths will be resolved progressively less clearly by the 
satellite sensors. This is obviously a pessimistic estimate of the vertical reso- 
lution as the results from 42 and 28 km appear to be quite independent. 

3. The Observed Temperature-Ozone Covariance 

The spectral analysis was performed by calculating the auto- and cross- 
correlation functions of the appropriate pairs of time series, applying a Parzen 
lag window, and then taking the Fourier transforms to give the auto- and cross- 
spectra (Jenldns and Watts, 1969). 

The occasional "wild" point was removed from the temperature and ozone 
data by rejecting all observations differing from the mean values by more than 
three sample standard deviations. The gaps, due to rejected observations and 
missing data, were filled by linear interpolation between values at either end of 
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the gap. The amount of missing data, including "wild" points, was 17% for 
temperatures and 31% for ozone. 

The temperature data is derived from a global estimate for one day, Green- 
wich time, while the ozone data is obtained from individual orbits and there- 
fore corresponds to local noon at 172 E. There is, therefore, an effective time 
shift of approximately 12 hours between the temperature and ozone time series. 
An attempt was made to remove this shift by linear interpolation between suc- 
cessive temperature observations. Any changes in the spectra were masked by 
the large fluctuations introduced in the coherence and phase spectra at the 
higher frequencies. This presumably occurred because the linear interpolation 
is equivalent to a two -point equally weighted running mean with its attendant 
sideband problems in the frequency domain. Consequently, the data has not been 
corrected for this 12 hour shift. It causes a phase displacement which is linear 
with frequency, 0 ° at 0 cycles per day (cpd) and -90® at 0. 5 cpd. At the frequen- 
cies of interest In this study, (0.2 cpd and less), the phase shift is comparable 
with the phase confidence limits. Much of the subsequent discussion is based 
on the limiting value of phase as the frequency tends to zero, so the phase shift 
will not be significant. 

The correlation functions were calculated for a maximum lag of 20 giving 
a maximum resolvable period of 40 days - a minimum resolvable frequency of 
0. 025 cpd. 

The 1971 austral spring transition is included in the data selected for anal- 
ysis. During this period the data will not be statistically stationary so that some 
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appropriate criterion was needed to distinguish spring from summer. Fortunately, 
an analysis of southern hemisphere temperature waves based on Nimbus 4 data 
was available (Harwood, 1975). Harwood found that by mid-October, 1971, there 
was a reduction in the amplitude of wavenumber 2 at 2 mb and a cessation of its 
eastward phase progression. The phase velocity of wavenumber 1 was eastward 
before mid-October, and westward after that time. Reversals in the meridional 
temperature gradient occurred at 2 mb about 28 September and at 10 mb about 
23 October. Accordingly 25 October was selected as a date by which the summer 
circulation was established in the stratosphere. The summer data is thus likely 
to be stationary although the spring data will not be. February 29, 1972, was 
selected as the end of summer. The division into spring and summer was made 
after the gap-filling procedure was completed. 

4. Planetary-Wave Perturbation Temperature and Velocity 

The observed phase spectrum (see section 6) often appears to have a fairly 
simple behavior - a phase of 0° over a wide range of frequencies, or 180° at low 
frequencies with a drift to 0° as the frequency increases. 

A simple model has been derived to see whether a comparable phase spec- 
trum is produced by a planetary wave mode with a given zonal wave number. 

A simple harmonic planetary wave is assumed to propagate vertically upwards 
through an atmosphere whose temperature and wind speed do not change with 
altitude. The ozone content is defined by the mass mixing ratio. The mean 
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mixing ratio is assumed to vary in the vertical and north-south directions (bai 
not in the east-west direction) to allow for vertical and meridional advection. 

Tlie finite temperature dependence of the photochemical equilibrium ozone 
concentration at 2 mb (Barnett, Houghton and Pyle, 1975) also requires some 
allowance for meridional temperature advection. A finite meridional tempera- 
ture gradient is therefore included, although the mean zonal wind is assumed to 
be constant at all heights. This blatant contradiction of the thermal wind equa- 
tion is necessary to avoid using Liouville-Green ("WKB") or numerical methods 
to cope with the variation of wave phase speed in the (vertical) direction of 
propagation which would result from a vertical variation of wind speed. 

The physical variables are represented by the following symbols: 

A j = amplitude of the perturbation stream function 
B = temperature dependence of equilibrium ozone concentration 
c = zonal phase velocity of the wave 
fg = Coriolis parameter for 45 S 
II = scale height 

k,l,m = eastward, northward and upward wave numbers (radians m‘^) 

N = buoyancy frequency 
P = atmospheric constant (2HN^/f^j) 

Q = ratio of horizontal to vertical advection terms 
R = gas constant for air 
R,. = radius of the earth 
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Tq = mean temperature 
T' = wave perturbation temperature 
Tj - total perturbation temperature 
u Q = mean zonal wind speed 

u' ,v ',w' = eastward, northward and upward wave perturbation velocity components 
x,y,z = eastward, northward and upward displacement coordinates 
X = ozone mass mixing ratio 

X^,X',Xj = mean, perturbation and total pex'turbation mass mixing ratios 

P ~ latitude dependence of Coriolis parameter 
3 

^ 3z 

7} = temperature coefficient of equilibrium ozone mass mixing ratio 

(j> = cross “Spectral phase 

\p' = wave perturbation stream function 

u) = wave angr.lur frequency 

In mathematical form the assumptions of the model are: 
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The wave properties can be deduced by substituting the above assumptions 
in the equation for the quasi-geostrophic perturbation potential vorticlty on a 
mid-latitude ^3-plane in the absence of heating and the equation for the thermo- 
dynamic balance (Charney 1973, Holton 1975). 



We assume the perturbation stream function appropriate to a simple 
harmonic plane wave 

4 /' = Aj exp [i(kx+£y + mz — cot) + z/2H] (3) 


It is convenient to define 


T = im + 1/2H 


so that 




Substituting (3) and (4) in (1) we obtain 


|(o;-kUo) [(k^ +K2) + (f2/N2)(m2 + 1/4H2)] -i-|3k| = 0 


( 4 ) 


( 5 ) 
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Eq. (5) has a non- trivial solution if 


w = kiig — j3k/[(k^ + 2^) + (fg/N^) (m^ + 1/4H^)] (6) 

Eq. (6) is the required dispersion equation from which we derive the horizontal 
^hase velocity 


C= w/k = Ug -|3/[(k2 + + i/4h2)] (7) 

The vertical wave number is given by 

m2 = (N2/f2 ) _ c) _ (k2 + g2 )] _ 1 /4 h2 (8) 

The vertical group velocity is 


|^= 2j3r2N^mk/[N^(k2 +22) + f2(,n2 + i/4h2)] 


(9) 


The vertical group velocity gives the appropriate sign for the vertical 
wave number (m) in (8). It will be assumed that energy is propagating upwards 
through the region; the vertical (m) and zonal (k) wave numbers must then have 
the same sign. 

From die stream function we obtain the horizontal perturbation velocity 
components and the perturbation temperature: 


u' = -x^' 
v' = ikxp' 


( 10 ) 

( 11 ) 
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From the thermodynamic equation (2) and from (4) we have the vertical 


perturbation velocity 

w'= i(7f^/N^)(6j-kuQ)i//' 

= i(R/N2)(o;-kUo)T' 

Eqs. (4), (11), and (12) give the relationship between the meridional per- 
turbation velocity and the perturbation temperature 


(13) 

(14) 


V = i(kR/7fo)T. 


(15) 


5. A Simple Model For Temperatrre-Ozone Covariance 
The rate of change of ozone mass mixing ratio is 


dX _ 3X 


_ OA . 1 9X . I 

— + V -T + W 


dt 


at 


ay 


ax 

az 


(16) 


Substituting (14) and (15) in (16) we have 




(17) 


Eq. (17) contains the basic principle of this analysis — that the phase difference 
between the ozone and temperature variations depends on a combination of atmo- 
spheric and wave parameters. It should therefore be possible to deduce some 
atmospheric parameters from the frequency dependence of the phase difference 



between temperature and ozone perturbations, supported by assumptions or 
measurements of the wave parameters. 

The time-dependent term in (16) and (17) would be zero if the ozone mass 
'mixing ratio was a conservative property of the atmosphere. Even if photo- 
chemical equilibrium is assumed there still remains the temperature dependence 
of the equilibrium mixing ratio. From the discussion by Barnett, Houghton and 
Pyle (1975) on the temperature dependence of the ozone concentration near the 
stratopause it is estimated that near 2 mb (42 Ian) the mixing ratio will behave 
like exp (B/T), where B ~1000 K. Consequently, the temperature coefficient of 
the ozone mass mixing ratio is 


1 ax 

^ “ X DT 


= ~B/T2 


Hence 


Let the mixing ratio be 


ax ^ ax dT 
at 3T dt 



X = Xg + Xj 


(18) 


(19) 


where 


Xj ~ exp (— icot) 


(20) 
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Therefore 


(21) 

Combining (11), (13), (16), (19), (20) and (21), and neglecting the spatial gradient 
of the perturbation mixing ratio, we have 

3X/, 

The time-dependent temperature term will contain a direct contribution from 
the wave perturbation temperature and another from the advective heating. The 
temperature is 

T = Tq + Tj 


where 


Therefore 


Tj ~ exp (— icjt) 


From (12) 


dT 9T' , 

"dT^ir 


ar 

at 


“i^^(7fo/R) 0' 


(23) 
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Consequently 


dt 


— icoTj 


= -ico(TfQ/R)i//' + i-^ k\p‘ 


and 


Tj = (7fo/R)\!''-(k/co)-^f 

\ 

Substituting (24) in (22) we obtain 


(24) 


(25) 


-itoXi = 


3T„ 


-iw7(fp/R)i//' + ki//' 


+ i-^k^' + i7(fo/N2){w-kUo)i//' 


so that 


X, = [7X^+(X^ + X,,)/co + 7\(Uo-c)/co] t//' 


where 


(26) 


- 7?Xo(fo/R) 
3X„ 


X,. = -k- 


Dy 


X. = -k77X 


9T. 


0 0y 


X 


V 



3z 
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These four terms represent the changes in ozone mass mixing ratio due to, 
respectively, the direct temperature change from wave perturbation, horizontal 
advection of the mixing ratio gradient, the temperature change from horizontal 
temperature advection, and vertical advection of the mixing ratio gradient 
Similarly the temperature change (25) is the sum of two terms 

T, = + (27) 

where T^ = fg/E, a direct perturbation term and Tj^ = -k 9Tg/3y, an advec*- 
tive term. 

The covariance of temperature and mixing ratio is then 
C(T,X.co) = TjXf 

(28) 

= [tT^ + + (X^ + Xj,)/w + (Ug - c)/co] 

The perturbations, Tj and Xj , are simple harmonic terms so that the coherence 
spectrum is of little interest, being perfect (1, 0) at all frequencies. However the 
phase spectrum (^(co)) does contain information. 

0(cu) = tan“i (C, (T. X, oj)/C^ (T, X, w) ) (29) 

where 

= real part of C(T,X,co) 

Cj = imaginary part of C(T, X,cj). 
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In a homogeneous atmosphere 


SX, 


3z 


0 


and 


9y 


= 0 


and tlie numerator of (28) is zero. The sign of tlie denominator depends only on 
the sign of tj, which is a negative number. Consequently, 0 (cd)= 18O°. The 
correlation coefficient between the two simple harmonic variables is then 


cos (0(to) ) = ~1, 0, 

Wlierever 90°, the correlation coefficient will be very small, even if the 
coherence is perfect. Consequently, it is inadvisable to use the correlation 
coefficient to express the relationsMp between quasi-periodic variables; an 
insignlfi,cant correlation coefficient does not necessarily imply a lack of 
dependence — it may merely indicate a phase difference close to 90°. 


6. Comparison of the Model and Observed Covariances 

In relating the model to observed covariances, the number of free param- 
eters has been minimized by taJdng the mean temperature and the mean mass 
mixing ratio and their respective gradients from the data itself. The values 
used are given in Table 1. Wliere magnitudes did not differ greatly from each 
other, a representative value was selected to simplify comparsions. For 
instance, the summer meridional temperature gradients at 2 and 20 mb were 
-2. 5 and -i-3. 0 K over 10° of latitude, while the 20 mb spring gradient was +1. 3K 
over 10° of latitude. A magnitude of 2 K over 10° of latitude was therefore 

chosen, witli the appropriate signs, for each of these values. 
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The 2 mb observations represent conditions of approximate photochemical 
equilibritim, with a linown temperature dependence of tlie ozone concentration | 
(Barnett, Houghton and Pyle, 1975). Conversely, the equilibrium time at 20 mb 
is of the order of months (Nicolet, 1975) so that on the time scales relevant 
here, only advective changes need be considered. Tliis is achieved by setting to 
zero the value of the temperature coefficient of the ozone mass mixing ratio at 
20 mb. : 

The selection of suitable zonal winds is also difficult as there is consider- 
able variation in tlie spring, and some variation during the four summer months. 
Calculations with a wide range of model parameters contirmed that the dominant 
effect of the wind is to determine when a given wave mode becomes evanescent, 
through the term in (8). The wind speed also appears in the (Uq-c) term 

of (28). As tlie wind speed is always associated with wave number and with the 
independent variable a> in 

(Uq - c) = (Ug - m/k), 

it is possible to qualitatively estimate the effect of a wind speed change from 
the consequent changes in wand k necessary to maintain a given value of (Uq-c). 
The model assumes a typical wind speed of 20 m sec'^ , the direction being 
eastward in the spring and westward in the summer at botli 2 and 20 mb. 

In this preliminary study, no observational estimate of the meridional wave 
number was sought, so a value of zero was selected, making this model 
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consistent with that of Hirota (1971), Holton (1975) discusses the problems of 
meridional scale and polar boundary values in the (3 -plane approximation. 

The observed 2 mb spring phase spectrum (Fig. 1) tends to a limiting value 
of about 215® at low frequencies. This behavior is consistent with the model and 
Suggests that the limiting phase angle (as the frequency decreases to zero) may 
have some significance, despite the fact that the periods greater than 40 days 
are not directly observable. From (29) it can be shown that 


where 


Limit <^(a)) = tan'^ 2mnH/( (PQ/iu) - 1) 


(30) 


m„ = Limit m (equation 8) 

UJ-J-O 

P = 2HN2/fg 



Tlie product PQ is independent of the wave number and comprises a constant P 
(“6 X 10“^ m sec'^ in this model) and a dimensionless number (Q) wliich is die 
ratio of die horizontal to vertical advection terms. 

The limiting vertical wave number (m^ ) does depend on the zonal wave 
number* through (8). However, the consequent variation in phase is not great, 
as shown in Table 2. 
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Eq. (30) sliows that the limiting phase is controlled by the ratio (Q) of 
horizontal to vertical advection terms. For instance, if there is no vertical 
ozone gradient, the limiting phase angle will be 0 or 180 , depending on the 
reLative magnitudes of the horizontal temperature and ozone advection terms. 

Tills is in agreement with the phase difference between meridional displacement 
and perturbation temperature calculated from (15). Around 20 mb and below, Q 
“will simply be the ratio of horizontal to vertical ozone gradients as the equilib- 
rium time is so long. 

The summer phase spectrum at 2 mb (Fig. 2) is not signilicantly different 
from 180° for frequencies less than 0.15 cpd. Thereafter, the phase decreases, 
but the phase fluctuations increase as tlie coherence decreases. This behavior 
is also consistent with the model, although the 180° phase angle there corresponds 
to the evanescent model. 

Hirota (1975) reported a westward moving vertically propagating wave 
number 1 mode in the southern hemisphere for summer 1972-3, Tliis is con- 
sistent with the trend in Harwood’s (1975) analysis for October 1971. Hirota 
also reported eastward moving waves with wave numbers 2 and 3 and suggested 
that these might be evanescent modes as (u^-c) was negative. The pliase speed 
of the westward mode was -28 m sec'^ at 30 S and -16 m see‘‘ at 60 S. It 
woxtld therefore be about -22 m sec'* at 44 S but the mean zonal wind is not 
given, so no comparison can be made with this model, but if the wind were com- 
parable with that assumed for this model, the wave mode would be close to evanes- 
cence (and therefore close to the maximum response to a given amplitude of lower 
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boundary forcing — see Hlrota, 1971), However, the experimentally observed 
low frequency phase is best explained by an evanescent mode in the model as the 
propagating modes have phase differences of about 90“ near their low frequency 
limit. 

The 20 mb spring phase spectrum (Fig. 3) Is not significantly different at 
low frequencies from that predicted by the model. The 20 mb summer phase 
Spectrum is not considered, as the coherence spectrum nowhere reached the 95;% 
significance level, presumably because the westward winds suppressed propaga- 
tion. 

7. Conclusion 

In principle, then, ozone measurements by balloon from one location could 
provide an estimate of the meridional ozone gradient, provided there are enough 
temper'ature observations to determine tlie meridional temperature gradient. 

The feasibility of this approach will need more extensive analysis of (30) with 
various gradients and due allowance for the confidence limits of observable 
phase. The latter depend in turn on the coherence spectrum. The coherence and 
phase are amplitude-independent and give no Indication of the amount of energy 
at the various frequencies; the cross-amplitude spectrum is therefore essential 
to any discussion of relative energy contributions of the various wave modes. 

The model will then require modification to include the dependence of forced 
amplitude on the vertical wave number (Hirota, 1971) and consequent assmnptions 
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about the spatial and temporal spectra of the forcing perturbation at the lower 
boundary. i 

It is also necessary to confirm the assumption that the temperature and 
ozone perturbations are due to the same planetary wave. Tliis study must await 
the availability of ozone data for several seasons around a latitude circle. 

The summer phase spectra and Hirota’s (1975) results imply that evanes- 
cent modes need to be considered more closely; this will require a more care- 

I 

ful examination of boundary conditions and the direction of energy propagation 
than has been assumed here for propagating modes. 

A further and highly desirable extension of this study is to see whether 
the principle exemplified by (30) can be extended to total ozone measurements 
made by a ground-based u. v. spectrophotometer at one location. 
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Table 1 


Observed Mean Values Used as Model Parameters 


Heiglit 

1 

1 

Season 

3y 

ay 

9Xq 

dz 

42 kin (2 mb) 

spring 

-6 K/10° lat. 

+10-12 ,„-i 

-10-5 m-i 

42 km (2 mb) 

summer 

-2 

-10-13 

-10-5 

28 km (20 mb) 

spring 

+2 

+2 X 10-12 

-3 X 10-1° 

28 km (20 mb) 

summer 

-2 

+10-13 

-3 X 10-10 
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Table 2 


Low Frequency Limit Of Model Phase Spectrum in Spring 


Height 

Zonal Wave Number 

-3 

-2 

-1 

-1 

+2 

+3 

42 km (2 mb) 

157“ 

O 

CO 

144“ 

216° 

212“ 

203“ 

28 km (20 mb) 

356“ 

354“ 

353“ 

7“ 

6“ 

4“ 


Positive zonal wave numbers indicate eastward-moving waves; negative zonal wave 
numbers indicate westward-moving waves. 
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Figure la Data for 2 mb, spring 1971. Power spectra for temperature (T) and ozone mass 
mixing ratio (0 ), and the coherence spectrum (C). Power is in decibels rela- 
tive to the total variance about the mean. The 95 ^ significance level is for 
cohe rence. 
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Figure lb Data for 2 mb, spring 1971. The obsei'ved phase spectrum is marked X, 
some 95*? confidence limits are indicated. The curves labelled +n (or -n) 
are the spectra calculated from the model for eastward (or westward) 
travelling modes with wave number n. Evanescence is indicated by the 
dashed line. 
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Figure 2a Data for 2mb, summer 1971-2. Power spectra for temperature (T) and ozone 
mass mixing ratio (02>> the coherence spectrum (C). Powder is in decibels 
relative to the total variance alxmt the mean. The 95% significance level is 
for coherence. 



360 



X- f -3 

•o «- S 

« o 2 

1 i ”2 .S 

2 — js -a 

—• fj o 

c 9 to J2 
S "S OJ to 
? ^ «s rt 

U ^ -Q 

Ip °i 

“ I '2 >> 

0) o •“ 

rft — 


$ -a ^ "O 

tjas 

1 2 -3 .2 

o “o 5 lu 

c e u 

2 £ c 

H 2 ^ S 

«i - S 

<m‘ 2 £ £ 

' •? o 5 

.5 


Q W C 

0) c i2 u 

S| 3i^ 

§ "S rt £ 

^ i3 2 
£ in o 5 

«S t ^ ^ 

j= 

o S o> .t: 

CO Is ^ 

■3-0 0^ 

5 S 5 


29 



00Z+ 


33N3U3H03 





IBP) UiMOd 


30 


RhTKODUClBll.in (Ji 
t)liI«LSAL PAGE 16 Pi^ 


Hgure 3a Data for 20 mb, spring 1971. Power spectra for temperature (Ti and ozone mass 
mixing ratio (0 ), and the coherence spectrum (C). Power is in decibels relative 
to the total variance about the mean. The 95% significance level is for coherence. 




